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Abstract 

We assess local-time variations of epithermal-neutron count rates measured by the Lunar Prospector Neutron Spectrometer. We 
investigate the nature of these variations and find no evidence to support the idea that such variations are caused by diurnal varia¬ 
tions of hydrogen concentration across the lunar surface. Rather we find an anticorrelation between instrumental temperature and 
epithermal-neutron count rate. We have also found that the measured counts are dependent on the temperatures of the top decime¬ 
ters of the lunar subsurface as constrained by the Lunar Reconnaissance Orbiter Diviner Lunar Radiometer Experiment temperature 
measurements. Finally, we have made the first measurement of the effective leakage depth for epithermal-neutrons of ~20 cm. 
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1. Introduction 


The mapping of hydrogen on the near surface of the Moon 
has been one of the focal points of the exploration of the Moon. 
Neutron spectroscopy has been central to this e ndeavor and al 
lows one to probe the top few decimeters (e.g. [Feldman et al 


1998t [Lawrence et al.j [2006[ ). Neutron spectrometers detect 


leakage neutrons resulting from the interaction of galactic cos¬ 
mic rays (GCR) and the chemical elements present in the top 
layers of the planetary subsurface. These neutrons lose energy 
by interacting with the surrounding lunar material, and different 
elastic and inelastic scattering cross-sections produce different 
amounts of moderation and capture. Since the Moon has no 
atmosphere, cosmic rays reach the surface without being ab¬ 
sorbed and neutrons can leak from the near subsurface without 
being scattered further or absorbed. In addition, the Moon has a 
characteristic composition pattern (dichotomy between Fe and 
Ti-rich maria and Al-rich highlands), to which neutron mea¬ 
surements are highly sensitive. Finally, from returned samples, 
the lunar regolith is recog nized to be hydrogen poor at equa¬ 
torial latitudes (<100 ppm,[Heiken et al.[[1991[). This a priori 
information has played a rnajor role m me interpretation and 
understanding of the orbital neutron data. 

According to their kinetic energy, £’kin, neutrons are typi¬ 
cally divided into three different energy bands: fast, epither¬ 
mal and thermal. These correspond to £’kin > 0.5 MeV, 0.4 eV 
< £'kin < 0.5 MeV and < 0.4 eV respectively. Briefiy, 
thermal neutrons are particularly sensitive to Fe, Ti, Gd and Sm 
in the absence of H ( [Elphic et al.[ [1998[ [2000[). Epithermal- 
neutrons are strongly moderated by H ([Feldman et al.| [1993J 
and therefore provide a robust measure of H concentrations 
and their spatial variations. Multiple studies have reported the 
presence and concentrations of H enhancements at the lunar 
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poles dFeldm an et al.[ [1998[ [Eke et al.[ [20091 [Teodoro et al. 
2010[ 2014D as well as^ co ncentrations within non-polar re- 
gions ( [Lawrence et al.[[MT5] ). 

In this article we study the daily variations of epithermal- 
neutron count rates from the Lunar Prospector neutron spec¬ 
trometer (LPNS) to inquire if such potential variations are sen¬ 
sitive to daily variations in H concentrations. This study is 
prompted, in part, by a report of time-variable concentrations 
of surficial (top tens of micrometers) H 2 O/OH that have been 
observed with spectral reflectance data ([Sunshine et al.[[2009[). 
If such time variable H conc entrations extend to m acroscopic 
depths of multiple mm to cm ([Lawrence et al.[|2011[), then such 
variations could in principle be observed witn orbital neutron 
data. Preliminary analysis of data from the Lunar Exploration 
Neutron Detector (LEND) on the NASA Lunar Reconnaissance 
Orbiter (LRO) mission has shown a diurnal count rate variation, 
which has been interpre t ed as daily variations in H concentra¬ 
tions (ILivengood et al.[ [2014J. Such findings, if confirmed, 
could lead to important consequences to our understanding of 
the hydrogen distribution at the lunar surface. Here, we will 
show that LPNS epithermal-neutron count rates also show vari¬ 
ations with solar local time, but that these variations are anti 
correlated with variations in instrumental temperature. Eur- 
thermore, the LPNS daily variations show a dependence on the 
temperatures at the top few decimeters of the lunar subsur face 
as predicted by the nu merical work of [Little et al.[([2003[) and 
Lawrence et al.[([2006[). These findings suggest that rather than 
reflecting diurnal changes in the hydrogen content of the lunar 
regolith, the temporal fluctuations in the count rates are due to 
small residual systematic effects in the data reduction. 


In Section [^we present the datasets used throughout the 
article while infection we perform the analysis of the ep¬ 
ithermal dataset. We finish the article in Section [4] with our 
conclusions and discussion. 
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2. Datasets 


In this section we will briefly detail the main characteristics 
of the epithermal-neutron dataset gathered by the LPNS, which 
is in the public domain at the Planetary Data System (PDS) 
(http : //pds-geo sciences . wust 1. edu/missions/' 
lunarp/ ). We will also use temperature measurements of the 
lunar surface collected by the Diviner Lunar Radiometer on 
Lunar Reconnaissance Orbiter. These measurements are also 
available at PDS (http ://pds-qeosciences . wust 1. 
edu/missions/ Iro/diviner . htm). 

2.1. LPNS Counting data 

The Lunar Prospector Neutron Spectrometer, designed and 
built at the Los Alamos National Laboratory, consists of two gas 
proportional counters containing helium-3. The neutrons that 
are captured by helium-3 nuclei within the detectors produce 
a unique energy signature that is then detected and counted. 
The two sensors have different covers: one is covered in cad¬ 
mium, the other in tin. The cadmium cover screens out thermal 
neutrons while the tin counterpart does not. Differences in the 
counts between the two tubes thus quantify the detected num¬ 
ber of thermal neutrons. The neutron spectrometer was located 
at the end of a 2.5 m boom, thus reducing the background from 
the main body of the spacec raft. For further detai ls regarding 
the neutron spectrometer see|Feldman et al.|(|20041). 

The epithermal-neutron counting rate and spacecraft alti¬ 
tude time series are shown in Figure along with other pa¬ 
rameters (neutron sensor temperature and pulse height peak- 
channel position) used in this study. All data shown in Figured 
are coloured based on the differe nt mission phases describea 
in Table 2 of|Maurice et al.|(l2004|). Green shows data from the 
High 1 phase (average altitude ~1U0 km), cyan shows data from 
the High 2 phase (average altitude -100 km, flipped spacecraft 
spin axis), magenta shows data from the Low 1 phase (average 
altitude -40 km), and violet shows data from the Low 2 phase 
(average altitude -30 km). The origin of time was set at 1st of 
January 1998 

Hereafter we will use data within the latitude range [-55°, 
55°], coming either from the High 1 (hereafter High Altitude) 
or Low 2 (hereafter. Low-altitude) phases, respectively. A third 
dataset with all the measurements within the same latitude range 
will be termed Overall. The number of 8s measurements within 
Overall, High altitude and Low-altitude datasets are 2,851,632, 
1,324,569 and 939,668, respectively. With this choice of lati¬ 
tudes we are considering those parts of the lunar surface whose 
epithermal signal is not dominated by the hydrogen enhance¬ 
ments in the top few decimeters of regolith. The individual 
count rates were normalized to altitude 30 km dMaurice et al.l 
[2^ . - 

2.2. Lunar Prospector Instrumental Temperature and Energy 
spectrum information 

The LPNS included a temperature sensor. Its time-series 
variation is presented in panel c) of Figure As with the count 
rate data we will use data within |latitude| <55°. 

The LPNS returned pulse height spectra of the ^He neutron 
captures, and a k^ feature is the 764 keV full-energy capture 
peak (see Figure [^. Panel d) of Figure shows the Cd ^He 
sensor peak-channel position versus time. The large jump in 
peak-channel position around time -300 days coincides with 
an increase in the detector’s high voltage. 
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Figure 1: a) Measured 8s counting rates and b) spacecraft altitude versus time. 
We are only considering data in which |latitude| < 55 . The colour scheme 
denotes different phases of the Lunar Prospector mission: : High 1 (green), 
High 2 (cyan), magenta (Low 1), violet (Low 2). The origin of the global time 
was set at 1st of Januar y 1998. The individua l counting rates were normalized 
to an altitude of 30 km ([Maurice et al.||2004J. c) Neutron sensor temperature, 
d) Peak-channel position ot the 704 keV neutron-capture peak in the LPNS 
pulse height spectra (see Figure]^. The large jump in peak channel-position 
around time ~300 days coincides with an increase in the sensor’s high voltage. 


2.3. Diviner/LRO temperature diurnal variation data 

The Diviner Lunar Radar Experiment on board LRO is an 
imaging multi-channel solar reflectance and infrared radiome¬ 
ter. This instrument maps the temperature of the lunar surface at 
< 500 meter horizontal scales. Diviner data sets are produced 
by the Diviner Science Team at the University of California, 
Los Angeles. 

Here we use temperature information as function of the lu¬ 
nar local time at different depths of the lunar subsurface. We 
will employ the Diviner measurements and model predicti ons 
constrained by such measurements (see |Siegler et al.|[201 1[ for 
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Figure 2: Upper left Diviner Temperature measurements at different latitudes function of the lunar local time. From the left to the right, and from the top to the 
bottom we show subsurface temperatures drawn from Diviner thermal m odels at 1, 2, 3, 4, 5, 10, 20 and 30 cm deep, respectively. These temperatures have been 
constrained by the Diviner observations using the methodology detailed in |Siegler et al.H20lT) . 


details). Upper left panel in Figure shows the Diviner mea¬ 
surements and the remaining eight panels show the model’s pre¬ 
dictions at depths of 1, 2, 3, 4, 5, 10, 20 and 30 cm. 

As expected, subsurface temperatures have a strong depen¬ 
dence on the local time, depth and latitude. At the surface, 
temperature peaks at the local midday. At deeper depths, the 
temperature maximum shifts towards later times of the day, as 
expected for a diurnal temperature wave. 


3. Analysis 


To quantify the variability of the epithermal count rate, cr(x, 
t), we introduce the stochastic overcount rate variable, dcr. 


6cr(x, t) = 


cr(x, t) - crf(x) cr(x, t) 


crf(x) 


crf(x) 


- 1 , 


( 1 ) 


where crf(x) = crj(lon.,lat.) is the fiducial count rate map at 
a given local time instant. Thus, the overcount rate field quan¬ 
tifies the change in counting rate, cr(x, t) - cr/(x), relative to 
a fiducial state, crfix). Hereafter, we will present results ob¬ 
tained using a fiducial state defined as the mean epithermal 


counting rate in a two-hour bin centered around 6 pm = -f90 
degrees from the sub-solar longitude on the Moon. Figure 
presents the fiducial count rate maps for the overall, high alti¬ 
tude and low-altitude datasets. cr(x, t) and crfix) are pixellated 
in 5°x 5° bins: crixijJk) and crfixij), where the spatial coor¬ 
dinate X is partitioned in latitude i e {1, • • •, 36} and longitude 
j e {1, • • •, 72), while the two-hour local time bins are repre¬ 
sented by k G {1, • • •, 12). In Figurewe present the averaged 
overcount rate, {Scr(Xij, tk))ij, in twomour local time bins cen¬ 
tered around 4 


{6cr{xij,t]fii)ij — 


- \ dcr(x,tk)W(Xij), 

^pixels (^/:) • 


( 2 ) 


where (• • •)/y denotes the average over all the pixels at a given 
instant 4, the sum is over all the spatial pixels (/, j), W(Xij) is a 
spatial “mask” applied to the data as the means to exclude the 
polar regions: 


W(xtj) = 


1, at latitudes g [-55°,55°], 
0, everywhere else. 


( 3 ) 
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Figure 3: Fiducial states adopted in this article, a) Overall data, b) High Alti¬ 
tude. c) Low-Altitude. The data is pixellized in 5°x 5° bins. 



Figure 4: Upper panel Average overcount rate per local time bin. The fidu¬ 
cial map is the mean counting rate at a given location over all the time series. 
The random errors in the measurements are smaller than 2.1x10“^, 5.2 xl0“^, 
3.5x10“^ for the overall (black), high altitude (green) and low-altitude (violet) 
data, respectively. Lower panel Average NS temperature per local time bin. 
The random errors in the measurements are smaller than 3.4x10“^, 8.2 xl0“^, 
5.2x10“^ °C for the overall (black), high altitude (green) and low-altitude (vio¬ 
let) data, respectively. The shaded regions of the diagrams represent the night- 
side. The x-axis tickmarks coincide with the local time bin centers. 


and 


^pixelsC^^) — ^ 


(4) 


The three derived time series (overall, high and low-altitude) 
show significant peak-to-peak variations, ~ 0.0273, given that 
the random errors are smaller than 3.0 x 10“^. 

In principal, these epithermal-neutron local-time variations 
may be due to variations in bulk hydrogen conc entrations across 
the lun ar surface, which is an idea suggested by|Livengood et al. 
(2014|) to explain similar variations in LRO/LUND data. How- 
ever, a comparison of LPNS and LEND data argue against such 
an interpretation. Specifically, even though the magnitude of 
the LPNS and LEND local-time variations are similar (2.7% 
for LPNS versus 2% for LEND), their qualitative behavior is 
fundamentally different. The LEND data show a single max¬ 
ima (6 local time [LT]) and minima (15 LT), whereas the LPNS 
data show two maxima (7 LT and 15 LT) and minima (12 LT 
and 24 LT). If there were true local-time hydrogen variations 
that caused a 2 to 2.7% variation in epithermal-neutron count 
rates, the measured local-time behavior should be similar in 
both datasets. The fact that their local-time behavior is dif¬ 
ferent suggests that other non-hydrogen systematic variations 
may be responsible for the observed local-time variations. Be¬ 
cause the two sets of measurements were made with different 
instruments on different spacecraft, such non-hydrogen varia¬ 
tions need not be identical in both datasets. Here, we inves¬ 
tigate possible sources of non-hydrogen systematic variations 
in the LPNS data. If we can show that such systematic effects 
can account for all or most of the local-time variations, then we 
can conclude that varying hydrogen concentrations are not the 
cause of these variations. Below, we discuss two such system¬ 
atic effects: variations in LPNS sensor temperature and varia¬ 
tions in lunar surface and subsurface temperatures. 


3.1. Epithermal-neutron count rate dependence on LPNS sen¬ 
sor temperature 

We first investigate if there is evidence that LPNS sensor 
temperatures might be related to systematic variations in the 
epithermal-neutron count rate that have not yet been removed 
in the prior data reduction procedures ([Maurice et al.| |2004|). 
Specifically, we examine if the varying sensor temperatures are 
related to gain variations (or a multiplicative shift) in the neu¬ 
tron sensor pulse height spectra. If such gain variations exist, 
we then ask if such variations are related to the 2.7% count 
rate variations. We note that the variations in the high and low- 
altitude data are qualitatively similar, so unless otherwise noted, 
we hereafter only carry out analyses for the low-altitude LPNS 
data. 

Inspection of Pigures[^ and[^ shows there is a likely cor¬ 
relation between the sensor temperature and the gain. Eigurej^ 
(upper) shows this correlation explicitly, where the temperature 
versus peak-channel position for the low-altitude data is plotted. 
It is clear that the gain of the sensor shows a positive, but non¬ 
linear variation with sensor temperature. The dynamic range of 
the variation is about 0.5 to 0.7 channels, or a 2 to 3% variation. 

To investigate if such a gain change could cause a count rate 
variation of 1 to 2%, Eigure (lower) shows a summed pulse 
height spectrum from the Cd ^He sensor. The LPNS pulse- 
height spectra were measured with 32 channels, so with such 
data it is difficult to see if a 0.5 channel difference in gain 
could cause a counting rate shift. We therefore fit the spectra 
with a gaussian plus polynomial between channels 10 and 30 
and plotted the function using 300 channel values instead of 
32 (red dashed line). To simulate the effect of a temperature 
dependent gain-shift, we shifted the spectra using the largest 
possible peak-position ratio allowed by the data (25.0/24.3). 
Both the original and gain shifted spectra are shown in dashed 
lines (original in red; gain shifted in blue). When one sums 
the counts between channels 15 and 30 (which are the chan- 
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Figure 5: Upper panel Low-altitude Peak-channel Position versus instrumen¬ 
tal temperature. Lower panel. Pulse height vs Channel from the LP Cd^He 
sensor. Original and gain shifted spectra are denoted in dashed lines: original 
in red and gain shifted in blue (see text for details). 


nels used by |Maurice et al.| (|2004|) for the count rate change 
measurements), there is a U. /% counting rate change. While 
this count rate change is smaller than the measured change, it 
is in the same direction as the measured data, such that lower 
gain data (which corresponds to lower instrument temperature) 
shows a larger counting rate, and therefore demonstrates that 
a varying gain caused by a varying temperature can result in 
varying LPNS epithermal-neutron count rates. 


3.2. Neutron Spectrometer sensor temperature correction 

To remove the sensor temperature induced count rate varia¬ 
tion, we have implemented the following two-step procedure: i) 
find the linear correlation coefficient between low-altitude time- 
series count rates and their associated sensor temperatures, and 
ii) use the linear relationship between those two variables to 
“de-trend” the data. A least squares fit to the data produces 

cr = -0.0194 • Tsensor + 19.170, (5) 

where the sensor temperature, Tsensor, is expressed in °C and 
the measured count rate, cr, in This relationship is applied 
to produce the sensor temperature-corrected count rates, 

cr' =cr + 0.0194 • [T^ensor - T sensor) , (6) 


[-55’<Ial<55*] 




local time [hours] 


Figure 6: Low-altitude count rate versus local time. Upper panel Black de¬ 
notes the uncorrected count rate while violet and green represent sensor tem¬ 
perature corrected and sensor plus surface corrected, respectively. Lower panel 
Subsurface temperature corrected count rates. Black, orange, violet and green 
represent count rate data corrected with Diviner subsurface temperature models 
at 0, 5, 10 and 20 cm, respectively. The random errors in all the measurements 
are smaller than 3.4x10“^. 


where Tsensor denotes the mean sensor temperature. These cor¬ 
rected count rates lead to fiuctuations in the low-altitude over¬ 
count rates that are reduced from 2.73% to 1.60%, as shown in 
Table 3 and the violet line in the upper panel of Figure In 
the same panel of Figurej^we also present the uncorrectea data 
(green). 


3.3. Count rate correlation with lunar surface and subsurface 
temperature 

ILittle et al.| (|2003|) and |Lawrence et al.| (|2006]) used the 
Monte Carlo transport code MCNFX to study the transport of 
neutrons in typical lunar soils and found that thermal and ep¬ 
ithermal count rates are dependent on soil temperature. Here, 
we consider a variety of subsurface temperatures as functions 
of both lat itude and local tim e, constrained by Diviner mea¬ 
surements ( [Siegler et al. 2011 1. 


To take into account the effects of the lunar subsurface tem¬ 
perature we consider an “average soil” made up of equal frac¬ 
tions of FAN, Luna 20 (L20), Luna 24 (L24), ^ollo 11 (All) 
and Apollo 12 (A12). In panel a) of Figure [7] we reproduce 
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Soil type asou bsou 


FAN 

Lunar 20 (L20) 
Apollo 11 (All) 
Apollo 12 (A12) 
Lunar 24 (L24) 


20.80 0.00280 
20.70 0.00220 

20.40 0.00160 

20.08 0.00211 
19.97 0.00163 


Table 1: Various soil parameters, (asoil, bsoil), considered in this work. These 
parameters were obtained using a linear least squares fit to the points in panel a ) 
of FigurejT] For further details regarding these soils see |Lawrence et al.| ( |20061 . 


Depth 

^ depth 

^ depth 

0 cm 

19.46 

0.00036 

1 cm 

19.38 

0.00067 

2 cm 

19.19 

0.00143 

3 cm 

18.98 

0.00227 

4 cm 

18.98 

0.00227 

5 cm 

18.84 

0.00284 

10 cm 

18.88 

0.00266 

20 cm 

19.67 

-0.00056 

30 cm 

19.46 

0.00036 


Table 2: Equation ^linear least squares fit parameters, at the 

various depths useoto model the epithermal count rates versus soil temperature. 


some of the|Lawrence et al.|(|2006|) results for epithermal-neutron 
count rates trom these lunar soils. The dashed black line repre¬ 
sents an “average lunar soil”, with epithermal count rates hav¬ 
ing a temperature dependence given by 


cr = Usoii + bsoiiTsoii = 20.39 + 2.068 x 10 ^Tsou, (7) 

where the soil temperature Tsoii is expressed in Kelvin. In Table 
we present the various soil parameters, {usoiu bsoii)^ that are 
used in the calculation of asoU and bsou of Equation These 
parameters result from linear least squares fits to the points in 
panel a) of Figure [7] 

(Sub)surface temperatures are derived using the Diviner mea¬ 
surements as well as the models constrained by these measure¬ 
ments shown in Figure]^ To obtain the proper temperatures for 
an epithermal-neutron temperature correction, we interpolate in 
local time and latitude at a given depth. We fit a linear curve to 
the box-smoothed count rates versus (sub)surface temperature 
at various depths, such that 

depth - ^depth + ^depth'^soil, depth, ( 8 ) 

where Tgou, depth is expressed in Kelvin. The values of and 
^ depth various depths are shown in Table [21 In panel b) 

of Figure [7] we show the count rates corrected for instrumen¬ 
tal temperature, cr', versus (sub-)surface temperatures at vari¬ 
ous depths. The count rates are box-smoothed using the same 
number of observations per bin and are denoted by dotted lines. 
The corresponding solid lines represent the linear least squares 
fits at each depth. For the sake of the temperature corrections, 
the given depth represents the “effective” depth from which the 
neutrons are emitted and is an approximation to the actual sce¬ 
nario where detected neutrons are emitted from a distribution 
of depths. 


To match the “average lunar soil” count rate versus soil 
temperature predictions shown at the upper panel of Figure J7] 
and the LPNS-Diviner measurements of the same relationship 
shown at the lower panel of the same figure we introduce the 
following transform: 


affinal, depth ~ depth 

ip'soil ~ ^depth) ^ soil, depth ip depth ~ ^ soil) + 
ip soil ~ bdepth^ {T soil, depth ~ soil,depth^ • (9) 

Here, Tsou, depth is the average soil temperature at a given depth. 
Briefiy, the transform is made up of two parts: i) move the mea¬ 
sured point linear upwards is such a way that makes it coincide 
with the value of the model predictions at TsoU, depth (second and 
third terms of the rhs in the above equation), and ii) make the 
slope of the linear fit to the measured points match the model 
counterpart, which is denoted by the last term in Equation]^ 

In the upper panel of Figure we show the average over¬ 
count rate, (dcr), corrected for sensor temperature and surface 
temperature effects (black). In the lower panel of the same fig¬ 
ure we show the overcount rate as corrected for subsurface tem¬ 
peratures at different depths. Table 3 contains the peak-to-peak 
variation of the overcount rate for different subsurface temper¬ 
ature corrections, with the minimum variation occurring for the 
subsurface depth of 20 cm. This result shows that if we correct 
the LPNS epithermal-neutron count rates for temperature at an 
effective depth of 20 cm, we can minimize the local-time vari¬ 
ations. Because the remaining variation is quite small (<0.5%), 
we therefore find that there is no need to invoke local-time vary¬ 
ing hydrogen concentrations to explain most of the LPNS local¬ 
time epitherm al-neutron variations. We further note that|Little 
et al.| ( |2003| ) predicted that the “effective” depth for thermal- 
neutron emission is 30 g/cm^, which for a soil density of 1.5 
g/cm^, results in a physical depth of 20 cm. This is the same 
value as our data-constrained effective depth for epithermal- 
neutrons, and supports the interpretation that much of the mea¬ 
sured local-time neutron variations are caused by subsurface 
temperature variations. 


4. Discussion and Conclusions 


In the prior section, we have shown that the LPNS epithermal- 
neutron count rate data equatorward of 55° show a local-time 
variation that has a magnitude of 2.7%. We have found that cor¬ 
rections for varying LPNS sensor temperatures and lunar sub- 
surfac e temperatures can redu ce this 2.7% variation to less than 
0.5%. |Livengood et al.|(|2014|) have also found ~2% local-time 
variations m the LEND epithermal-neutron count rates. Their 
interpretation of this result was that these variations were cause 
by mobile hydrogen that deposits once the sun sets, since the 
surficial temperature is low enough to lead to a long residence 
time during the night. It then leaves the lunar surface at dawn 
as higher temperatures imply very short residence times (see 
ISchorghofer and Taylor||2007|) for a definition of residence time 
and Its dependence on temperature). Because we can quantita¬ 
tively account for all but 0.5% of the LPNS local-time variation, 
we find little-to-no support for the interpretation that varying 
hydrogen concentrations are the dominant cause of the local¬ 
time epithermal-neutron variations in either the LPNS or LEND 
datasets. 

Even though the remaining LPNS local-time variations do 
not exhibit the time dependence expected for diurnal variations 
suggested by Livengood et al. (2014), we can nevertheless in¬ 
vestigate the implications if the remaining 0.5% local-time vari¬ 
ations are due to time-variable hydrogen concentrations. Specif¬ 
ically, if a 0.5% diurnal change in epithermal-neutron counting 
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Figure 7: Panel a) This panel, as in|Lawrence et al.U2006), shows count rates 
versus surface temperatures. Each color denotes a ditterent soil: FAN (brown), 
Lunar 20 (L20, powder green), Apollo 11 (All, orange), Apollo (A12, violet) 
and Lunar 24 (L24, blue navy). The dashed black line denotes the “average lu¬ 
nar soil”. Panel b) The measured count rates versus surface and sub-subsurface 
temperature. Dotted lines represent the box-smoothed count rates. Solid lines 
denote the least squares fits to the smoothed data. For clarity the 5, 10 and 20 
cm curves are shifted upwards by 0.25, 0.50 and 0.75 s~^, respectively. 


LPNS 

sub-dataset 


{Scr)^ax 

Acr 

Overall 

-0.0160 

0.0075 

0.0235 

High altitude 

-0.0134 

0.0091 

0.0226 

Low-altitude 

-0.0144 

0.0130 

0.0273 

Low-altitude 




sensor temperature 




corrected 

-0.0075 

0.0075 

0.0150 

Low-altitude 




subsurface 




temperature 




corrected 




Depth 0 cm 

-0.0052 

0.0040 

0.0092 

Depth 1 cm 

-0.0044 

0.0048 

0.0092 

Depth 2 cm 

-0.0053 

0.0044 

0.0098 

Depth 3 cm 

-0.0056 

0.0052 

0.0108 

Depth 4 cm 

-0.0056 

0.0052 

0.0108 

Depth 5 cm 

-0.0042 

0.0057 

0.0098 

Depth 10 cm 

-0.0028 

0.0029 

0.0057 

Depth 20 cm 

-0.0027 

0.0022 

0.0049 

Depth 30 cm 

-0.0052 

0.0040 

0.0092 


Table 3: Peak-to-peak variation of the LPNS overcount rate variation, Acr = 
{Scr)min - {^cr)min, of the different datasets applied in this article. 


rate were the result of mobile water molecules, then how much 


water would this imply? We note that for this type of diurnally 
varying water, the water stratigraphy would likely take the form 
of a thin wet layer over a thicker dry layer. In this case, we use 
results from [Lawrence et al.| (|2011 ) who modeled epithermal- 
neutron count rates tor wet-over-dry layering and found that for 
top layers thinner than 5 to 10 cm, increasing amounts of hy¬ 
drogen would result in a epithermal-neutron count rate increase 
rather than decrease. Here, we therefore assume that hydrogen¬ 
bearing compounds are in an upper layer of limited thickness 
that participates in the diurnal condense /release cycle. We con¬ 
sider Figure 3 of [Lawrence et al.|(|201 lb and posit an upper, ex¬ 
changeable regolith layer thickness of 3 g/cm2. The epithermal 
count rate in LPNS would increase by ~0.5% if that upper layer 
contained ~1 wt% water-equivalent hydrogen, i.e., 0.03 g/cm^ 
of H 2 O, over dry regolith. This 30 mg/cm^ of exchangeable 
H 2 O corresponds to 0.03/18 moles/cm^ of H 2 O molecules, or 
10^^ molecules/cm^. If exchangeable and mobile, this number 
is the column density of exospheric water that would be present 
over the dayside. Assuming a rough dayside scale height of 100 
km (10^ cm), the near-surface density of H 2 O would be 10^"^ 
molecules/cm^. Preliminary reports show that observed levels 
of H 2 O and OH are at leas t 10 orders of magnitude below this 
value ([Benna et al.[ [2014J). Other mobility scenarios that do 
not involve the luna^r exosphere seem highly implausible. 

Alternatively, we suggest that the remaining 0.5% variation 
is likely to due other systematic variations that have not been 
fully taken into account. As one suggestion, we consider the 
fact that the subsurface temperature correction approximated a 
distribution of neutron emission depths using a single “effec¬ 
tive” neutron emission depth. While calculating the true dis¬ 
tribution of neutron emission depths is beyond the scope of 
this current study, we can carry out a slightly more compli¬ 
cated (and possibly realistic) subsurface temperature correction 
by assuming that the neutrons are uniformly emitted through¬ 
out the top 30 cm. With this assumption, we can reduce the 
prior 0.49% local-time variation to less than 0.25%. While this 
final result should be treated with some caution because a true 
depth-dependent emission function has not been derived, it does 
demonstrate that the remaining 0.49% variation shown in Ta¬ 
ble 3 can likely be reduced further without appealing to varying 
hydrogen concentrations. 

In summary, we have shown that the LPNS epithermal-neutron 
data do exhibit local-time variations that have a magnitude of 
2.7%. However, these variations can be understood as being 
due to systematic variations of the LPNS sensor temperature as 
well as lunar subsurface temperatures. An additional result of 
this work is that by showing the LPNS local-time variations can 
be minimized by a subsurface temperature at a depth of 20 cm, 
the LP neutron data have shown to be consistent with modeled 
subsurface temperatures, neutron-transport modeled tempera¬ 
ture dependencies, as well as the predicted effective depth of 
subsurface neutron emission. 
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